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ANALYSIS  OF  REGIONAL  BODYWAVE  PHASES 
FROM  EARTHQUAKES  IN  WESTERN  CHINA 


OBJECTIVE: 

The  purpose  of  this  study  is  to  improve  our  understanding  of  the  regional  crustal  waveguide 
phases  encompassed  by  P„  and  Pg  (often  termed  P.,),  by  modeling  regional  waveforms  in  eastern 
Asia.  Saikia  and  Burdick  (1991)  showed  that  deterministic  modeling  of  short-period  P^  can  provide 
a  good  fit  to  waveforms  at  ranges  of  200  to  420  km  from  NTS  explosions.  Zhao  and  Helmberger 
(1991)  have  demonstrated  similar  success  in  modeling  broad-band  Poj  from  the  Saguenay  earthquake 
recorded  at  HRV,  while  Burdick  et  al.  (1992)  modeled  regional  earthquake  recordings  at  the  IRIS 
station,  Garm. 

If  we  want  to  understand  the  development  and  propagation  of  P„,,  we  niust  have  observations 
at  a  number  of  ranges  from  the  source.  Unfortunately,  while  high-quality,  broad-band  seismometers 
are  now  available  in  China  and  the  former  Soviet  Union,  the  station  spacing  is  quite  sparse.  In  a 
previous  report,  we  presented  the  results  of  modeling  a  profile  of  earthquakes  recorded  at  the  Chinese 
Digital  Seismic  Station,  WMQ  (Barker,  1991;  Wu  and  Barker,  1992).  Since  the  earthquakes  have 
different  depths,  magnitudes  and  mechanisms,  this  is  not  the  reciprocal  problem  to  wave  propagation 
from  a  single  source  to  a  number  of  stations.  However,  by  modeling  simultaneously  waveforms 
from  each  of  these  earthquakes,  we  gain  an  understanding  not  only  of  the  propagation  of  specific 
crustal  phases  near  WMQ,  but  also  of  the  kinds  of  variation  observed  for  different  source  depths 
and  mechanisms. 

Since  the  conclusions  reached  in  modeling  waveforms  from  a  single  profile  of  earthquakes 
are  radiation-pattern  dependent,  we  now  present  the  results  of  similar  waveform  modeling  studies, 
but  using  profiles  of  earthquakes  along  three  different  azimuths  from  WMQ  (Figure  1;  Table  1). 
The  first  is  a  repeat  of  the  previous  study  of  earthquakes  southwest  of  WMQ.  These  events  occurred 
within  the  crust  in  the  Tienshan  region,  and  propagated  along  the  structural  trend  of  the  Tienshan 
and  the  Tarim  Basin.  A  second  profile  is  located  to  the  NNE  of  WMQ,  through  the  western  tip  of 
Mongolia  and  into  Siberia.  The  third  profile  is  SSE  of  WMQ,  including  earthquakes  in  the  vicinity 
of  the  Lop  Nor  Test  Site.  The  propagation  paths  cross  the  eastern  margin  of  the  Tienshan,  the  Tarim 
Basin  and,  at  the  largest  distances,  the  Altyn  Tagh  and  Qaidam  Basin. 


1 


RESEARCH  ACCOMPLISHED: 

SWPffifik 

As  a  firet  profile,  we  consider  earthquakes  along  a  line  SW  of  WMQ  (Figure  1  and  Table  1), 
from  the  Tarim  Basin  and  the  Tienshan  regions  of  western  China.  In  a  surface  wave  regionalization 
study  (Wu  and  Jones,  "Surface  wave  regionalization  and  tomography  in  China  and  its  vicinity”,  in 
Wu  and  Barker,  1992),  these  are  considered  to  be  within  the  same  structural  region,  so  lateral 
variations  in  crustal  structure  should  be  minimal.  These  are  shallow  crustal  earthquakes  (depths 
8-33  km)  with  thrust  mechanisms.  Broad-bnnd  seismograms  have  been  processed  to  facilitate 
comparison  with  synthetic  seismograms.  This  processing  includes  time  integration  (to  ground 
displacement)  and  a  high- pass  Butterworth  filter  (frequency  0.08  Hz)  to  reduce  low-frequency  drift 
in  the  synthetics.  In  this  study  we  are  concentrating  on  the  Pn-Pg  wavetrain,  so  only  the  vertical 
component  is  modeled. 

A  profile  of  the  vertical-component  waveforms  is  shown  on  the  left  side  of  Figure  2. 
Superimposed  on  the  waveforms  are  travel-time  curves  appropriate  for  various  P  and  S  phases  for 
a  source  at  30  km  depth  in  a  layered  velocity  structure  model  (discussed  below).  To  facilitate 
comparison  between  events,  the  waveforms  in  this  figure  and  the  next  have  been  band-pass  filtered 
from  0.5  -  2.0  Hz,  similar  to  the  WWSSN  short-period  band.  Also  time  shifts  have  been  applied 
to  three  of  the  records;  for  events  87005  (560  km)  and  87159  (1 175  km),  a  time  lead  of  2  sec  is 
used,  while  for  event  87 35 1  (422  km)  a  lag  of  3  sec  is  used.  These  may  reflect  errors  in  the  assumed 
origin  time  of  these  events,  or  simply  variations  due  to  source  depth.  The  first  50  sec  of  these 
waveforms,  which  includes  the  P„-Pg  wavetrain,  are  shown  on  the  left  side  of  Figure  3,  along  with 
travel-time  curves  for  selected  phases.  A  number  of  features  in  the  observed  waveforms  correlate 
with  some  of  these  predicted  arrivals.  In  particular,  for  the  closest  event,  P,  pP  and  S  may  be 
identified.  Beyond  400  km,  P„  and  P  may  be  identified,  but  the  Moho  reflection,  P^P,  is  not  a 
substantial  arrival.  In  fact,  for  these  mechanisms,  sP^  and  sP^P  may  be  seen  as  an  elongated  series 
of  arrivals  at  400-600  km  and  as  distinct  phases  at  1 175  km.  Many  other  arrivals  are  present  in  the 
observed  waveforms;  the  travel-time  curves  show  only  selected  arrivals  for  a  single  source  depth. 

Other  features  are  better  modeled  by  computing  synthetic  seismograms  for  the  appropriate 
range,  depth  and  mechanism  and  comparing  this  with  the  observed  waveform.  The  velocity  structure 
model  assumed  (Table  3,  Figure  4)  is  based  on  the  surface  wave  results  of  Feng  and  Teng  (1983), 
modified  so  that  the  travel-time  curves  provide  reasonable  agreement  to  observed  arrival  times  (as 
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in  Figure  3).  The  Moho  is  at  a  depth  of  56  km,  whiJe  a  mid-crustal  discontinuity  is  located  at  41 
km  depth.  In  the  figures  to  follow,  reflections  from  the  Moho  are  denoted  P^P,  while  those  from 
the  mid-crustal  discontinuity  are  denoted  PcP.  A  velocity  gradient  is  included  in  the  mantle  so  that 
P„  is  modeled  as  a  turning  ay  rather  than  as  a  head  wave.  Synthetic  seismograms  were  computed 
using  a  frequency- venumber  (F-K)  integration  technique  (Barker,  1984).  This  method  uses  the 
compound  matrix  modification  of  the  Haskell  layer  matrix  method  with  Filon  quadrature  over 
wavenu.nber.  Anelastic  attenuation  is  included  to  move  the  poles  off  of  the  real-/:  axis.  No 
'’•avenumber  filtering  is  imposed,  so  the  synthetics  include  S  waves  and  surface  waves  in  addition 
to  the  P  wavetrain.  These  are  computation-intensive  synthetics,  so  we  must  limit  the  frequency 
band  and  time  duration  (up  to  4  Hz,  5 12  sec  duration).  The  source  parameters  used  in  generating 
the  synthetics  are  listed  in  Table  2.  These  include  Harvard  CMT  mechanisms  (published  in  the 
PDF)  when  available;  otherwise  an  average  mechanism  is  assumed.  Source  comer  frequencies  and 
Butterworth  filter  parameters  are  chosen  to  give  the  best  agreement  between  data  and  synthetics. 
Source  depths  (again  from  the  PDF)  are  sometimes  questionable,  so  the  synthetics  are  computed 
at  10,  20  and  30  km  depths,  and  the  depth  closest  to  that  reported  for  an  event  is  used  in  the 
comparison.  Generalized  ray  theory  synthetics  (Helmberger  and  Harkrider,  1978)  were  used  to 
identify  specific  rays. 

Profiles  of  F-K  synthetics  for  a  source  depth  of  30  km  are  shown  on  the  right  sides  of  Figures 
2  and  3.  Although  some  wrap-around  is  apparent  at  the  beginning  of  the  traces,  P„  and  several  later 
arrivals  may  be  easily  identified.  The  synthetics  are  somewhat  simpler  than  the  observed  waveforms 
(compare  the  two  sides  of  Figure  3),  but  many  features  are  common.  For  example,  at  400  km  sP„ 
and  sPmP  interfere  to  generate  an  elongated  wavetrain  near  20  sec  reduced  time.  Although  not 
shown  in  the  travel-time  curves,  the  second  and  third  P-wave  reverberations  in  the  crust  also  arrive 
between  20-30  sec  (reduced  time)  at  this  range.  With  increasing  range,  sP^  becomes  the  dominant 
phase,  interfering  with  PmP  at  1175  km  range.  Higher-order  crustal  multiples  (PmPPmP.  SmPPmP^ 
etc.)  do  not  appear  to  play  a  dominant  role  in  either  the  observed  or  synthetic  waveforms  for  these 
earthquakes.  Certainly  the  strength  of  the  upgoing  S  wave  that  reflects  from  the  free  surface  is 
dependent  on  the  radiation  pattern,  and  in  this  profile  we  are  considering  only  eanhquakes  along  a 
single  azimuth  and  with  comparable  mechanisms.  For  near-surface  isotropic  sources  (explosions), 
we  would  expect  crustal  multiples  to  dominate  the  waveform  as  Burdick  el  al.  (1989)  found  for 
NTS.  This  is  an  example  of  how  radiation  pattern  can  cause  substantial  difference  in  the  generation 
of  the  high-frequency  P„-Pg  waveform,  and  may  be  exploited  as  a  discriminant. 
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Since  the  Pn-P,  waveforms  result  from  the  interference  of  a  number  of  phases  which  depart 
the  source  either  upward  or  downward,  it  is  instructive  to  see  how  this  interference  varies  with 
changes  in  source  depth.  Shown  in  Figures  5  -  9  are  observed  vertical-component  waveforms  for 
five  of  the  events  in  the  SW  profile,  along  with  F-K  synthetics  computed  for  10,  20  and  30  km 
source  depths.  With  the  exception  of  event  87279  (Figure  5),  the  traces  have  been  aligned  on  the 
P„  arrival  (87279  is  at  prc-citical  range,  so  is  aligned  on  P).  Upward  departing  phases  (such  as  sP^) 
move  out  in  time  with  increasing  source  depth,  while  downward  departing  phases  (such  as  P^P) 
remain  stationary  or  move  in.  The  arrival  times  of  important  phases,  determined  from  generalized 
ray  synthetics,  are  suggested  by  the  lines  on  the  figures.  Of  course,  since  we  have  assumed  a  layered 
structure,  the  change  of  arrival  time  with  depth  is  not  really  linear;  the  lines  are  used  to  help  the 
reader  visually  interpolate  arrivals  between  the  depths  we  have  chosen  for  compulation.  Clearly, 
for  different  depths  and  different  ranges,  different  phases  interfere  to  form  the  arrivals  observed  on 
the  venical-component  seismograms.  For  example,  for  event  87279  (82  km.  Figure  5),  P^P  is  a 
relatively  minor  phase,  but  pPmP  and  (PcP)2  (a  double  reverberation  in  the  upper  crust)  interfere  at 
30  km  depth  to  produce  a  single  large-amplitude  arrival,  which  corresponds  to  the  largest  arrival 
in  the  observed  seismogram.  The  published  depth  for  this  event  is  32  km,  which  is  consistent  with 
the  depth  inferred  from  the  synthetics  (slightly  deeper  than  the  30  km  synthetics,  denoted  by  the 
arrow  in  Figure  5). 

For  event  87351  (422  km.  Figure  6),  crustal  phases  are  well  separated  in  time,  resulting  in 
the  elongated  series  of  arrivals  observed  for  this  event.  If  the  depth  is  somewhat  greater  than  30 
km  (as  indicated),  arrivals  observed  at  about  1 8  sec  and  33  sec  may  be  interpreted  as  sP„  and  s(PmP)2, 
respectively.  The  large-amplitude,  late  arrival  in  the  synthetics  is  Sn  which,  as  usual,  is  substantially 
larger  in  the  synthetics  than  in  the  observed  waveform.  At  560  km  (event  87005,  Figure  7),  none 
of  the  computed  synthetics  matches  the  arrival  times  of  all  of  the  observed  phases,  but  from  the 
relative  moveout  of  P„,  P  and  sP„,  we  can  sec  that  a  source  depth  of  14-15  km  would  produce  an 
excellent  fit.  At  this  depth,  the  large  arrival  at  30  sec  results  from  sPeP  and  s(P„)2.  The  published 
depth  for  this  event  is  17  km.  On  the  other  hand,  forevent  98024a  (731  km.  Figure  8),  a  source 
depth  of  about  1 6  km  would  provide  a  better  fit  (particularly  for  sPeP  and  s(PmP)2)  than  the  published 
depth  of  30  km.  Finally,  for  event  87159  (1175  km.  Figure  9),  the  published  mechanism  is  clearly 
inconsistent  with  the  observed  P-wave  polarities  at  WMQ.  However,  since  the  crustal  phases  are 
well  separated  in  time  at  this  range,  we  interpret  that  the  source  must  have  been  shallower  than  the 
published  depth  of  10  km. 
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NNE  Profile 

Another  profile  of  eanhquakes  is  located  to  the  NNE  of  WMQ,  crossing  the  Altai  Mountains, 
through  the  western  tip  of  Mongolia  and  into  Siberia  (Figure  1;  Table  1).  In  this  region,  very  little 
has  been  published  on  the  crustal  structure.  Surface-wave  tomography  (Wu,  1993)  suggests  that 
crustal  structure  varies  slowly  along  this  profile.  Therefore,  we  have  developed  a  crustal  velocity 
structure  based  on  fitting  travel-time  curves  to  observed  seismograms.  Once  again,  venical-com- 
ponent  broad-band  seismograms  were  obtained  from  CSS,  integrated  to  ground  displacement,  and 
high-pass  filtered  at  0.08  Hz.  The  seismograms  recorded  at  WMQ  from  the  NNE  profile  of 
earthquakes  are  shown  on  the  left  side  of  Figure  10,  along  with  travel-time  curves  computed  for  a 
source  at  20  km  depth  in  the  velocity  structure  in  Table  2.  The  Pq-P,  portions  of  these  seismograms 
are  plotted  on  the  left  side  of  Figure  1 1 .  The  linear  raoveout  of  Pg  (Figure  1 1 )  suggests  a  very  low 
velocity  gradient  in  the  crust.  P„  goes  from  8.15  km/s  at  400-600  km  range  to  8.26  km/s  at  greater 
ranges,  constraining  the  mantle  gradient.  If  we  identify  several  observed  secondary  arrivals  with 
P^P,  the  velocity  at  the  base  of  the  crust  is  7.4  km/s.  Since  the  depths  of  some  of  the  earthquakes 
in  this  profile  are  uncertain,  and  we  have  little  independent  constraint  on  crustal  structure,  if  the 
source  depth  is  other  than  20  km  we  need  to  vary  the  crustal  thickness  in  order  to  preserve  the  time 
separation  between  P„  and  P,.  The  resulting  velocity  structure  models  are  plotted  in  Figure  12.  For 
sources  at  10  km  depth,  the  crustal  thickness  is  52  km,  with  a  slight  gradient  at  the  base  of  the  crust. 
For  20  km  depth,  the  crustal  thickness  is  59  km,  while  for  30  km  depth  we  assume  a  crustal  thickness 
of  62  km. 

Synthetic  seismograms  were  once  again  computed  by  the  frequency- wavenumber  integration 
technique,  for  the  source  parameters  listed  in  Table  3.  In  this  case,  the  mechanisms  are  predominantly 
strike  slip,  and  we  assume  an  increasing  normal-fault  component  to  the  northeast  (toward  Lake 
Baikal).  The  eanhquakes  are  small  enough  that  the  source  comer  frequency  has  no  effect  within 
the  band  computed,  so  the  synthetics  are  simply  the  impulse  response  computed  to4  Hz.  Butterworth 
low-pass  filters  (3.0  Hz)  are  applied  to  both  observed  and  synthetic  data.  The  profiles  of  synthetic 
seismograms  for  a  source  at  20  km  depth  are  plotted  on  the  right  sides  of  Figures  10  and  11.  For 
these  profiles,  the  observed  and  synthetic  seismograms  have  been  bandpass  filtered  (0.5  -  2.0  Hz) 
to  approximate  the  WWSSN  short-period  response.  Once  again,  although  the  synthetics  are  simpler 
than  the  observed  seismograms,  the  timing  and  character  of  P„,  P,  and  S,  arc  quite  adequately 
modeled  (Figure  10).  Surface  reflected  phases  provide  arrivals  between  P„  and  P,  (Figure  11). 
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The  comparisons  between  seismograms  observed  at  WMQ  for  each  event  and  synthetics 
computed  for  depths  of  10,  20  and  30  km  are  shown  in  Figures  13-17.  The  closest  event,  87279 
(82  km.  Figure  13),  is  the  same  event  included  at  the  closest  range  in  the  SW  profile.  In  this  case, 
however,  the  observed  and  synthetic  seismograms  have  different  filters,  and  the  synthetics  are 
computed  for  a  different  mechanism  and  different  structure  model.  For  this  mechanism,  the  largest 
arrival  at  about  16  sec  is  likely  the  S  wave.  The  next  largest  arrival  (at  20  sec)  is  pPmP,  suggesting 
that  the  source  depth  is,  once  again,  about  32  km  (denoted  by  the  arrow). 

Event  87261  is  at  a  range  of  388  km  and  event  88092  is  at  a  range  of  439  km.  Rather  than 
model  each  separately,  we  compare  each  with  synthetics  computed  for  a  range  of  400  km.  For 
event  87261  (Figure  14),  the  first  arrival  is  Pj,  P^P  is  small,  and  the  large  arrivals  at  about  15  sec 
begin  with  PcP  and  P  and  are  followed  by  a  series  of  crustal  multiples.  The  best  agreement  between 
observed  and  synthetic  waveforms  is  found  between  20  and  30  km  depth.  Comparing  the  same 
synthetics  with  the  observed  waveform  for  event  88092  (Figure  15),  we  see  that  P,  P^P  and  PcP 
are  very  emergent,  P  arrives  somewhat  earlier  than  forevent  87261 ,  and  the  train  of  crustal  multiples 
is  slightly  more  compact  in  time.  For  these  reasons,  we  interpret  the  depth  as  1 1  km,  which  is 
consistent  with  the  depth  of  10  km  listed  in  the  PDE. 

The  observed  waveform  for  event  88205  (586  km.  Figure  16)  is  rather  noisy  and  oscillatory. 
The  arrivals  in  the  synthetic  seismograms  are  much  more  distinct.  Nevenheless,  we  can  interpret 
the  large  arrival  at  about  22  sec  as  the  direct  P  wave,  the  smaller  arrival  preceding  it  as  PcP,  and 
the  small  first  arrival  as  P„.  With  some  imagination.  P^P  may  be  identified  and,  possibly,  pPmP- 
The  large  arrivals  following  the  direct  P  wave  are  most  compact  for  a  shallow  source,  and  become 
elongated  for  a  deeper  source.  These  pieces  of  evidence  suggest  that  the  source  of  this  event  is  at 
about  19  km  depth,  which  is  the  depth  listed  in  the  PDE.  The  Pn  portion  of  the  observed  waveform 
forevent  88182  (759  km,  Figure  17)  is  also  quite  noisy,  and  prevents  an  alignment  of  the  waveforms 
on  P„;  therefore  we  have  aligned  the  synthetics  on  the  P  wave.  pPmP  precedes  the  P  wave,  and  can 
be  seen  at  about  30  sec.  Following  P  are  the  crustal  multiples.  From  the  arrival  times  of  pPmP,  P 
and  p(PmP)2,  we  interpret  the  depth  as  16  km.  The  PDE  depth  was  listed  as  the  default  value  of  33 
km. 

SSE  Profile 

A  final  profile  of  earthquakes  recorded  at  WMQ  consists  of  events  located  to  the  SSE  (Figure 
1 ;  Table  1 ).  The  closest  event  (87279)  is  the  same  one  included  in  the  two  previous  profiles.  Three 


events  are  located  within  the  Tarim  Basin  in  the  vicinity  of  the  Lop  Nor  Test  Site.  Two  more  distant 
events  are  located  within  the  Altyn  Tagh  and  the  Qaidam  Basin.  Propagation  from  each  of  these 
events  crosses  the  Tarim  Basin  and  the  eastern  margin  of  the  Tienshan  before  arriving  at  WMQ. 
Vertical-component  broad-band  seismograms  were  obtained  from  CSS,  integrated  to  displacement, 
and  low-pass  filtered  (0.08  Hz).  For  the  profiles,  observed  and  synthetic  seismograms  have,  once 
again,  been  band-pass  filtered  (0.5  -  2.0  Hz).  The  seismograms  recorded  at  WMQ  from  the  SSE 
profile  of  earthquakes  are  shown  on  the  left  side  of  Figure  1 8,  along  with  travel-time  curves  computed 
for  a  source  at  20  km  depth  in  the  velocity  structure  in  Table  2.  The  P„-P,  portions  of  these  seis¬ 
mograms  are  plotted  on  the  left  side  of  Figure  19. 

The  velocity  structure  (Figure  20,  T  able  2)  is  modified  fit)m  models  of  surface  wave  dispersion 
in  the  Tarim  Basin  and  seismic  refraction  in  the  Qaidam  Basin  (SSB,  1986),  removing  thin  low- 
and  high-velocity  layers  and  adjusting  crustal  thickness  and  velocities  so  that  the  travel-time  curves 
are  in  agreement  with  the  observed  arrival  times.  As  noted  in  Table  1 ,  there  is  substantial  uncertainty 
in  the  depths  of  the  earthquakes  listed  in  the  PDE,  but  those  that  are  well  determined  fall  in  the 
21-32  km  depth  range.  The  travel-time  curves  are  plotted  assuming  that  most  of  the  earthquakes 
occurred  near  20  km  depth.  If  the  average  depth  of  the  sources  is  actually  deeper,  then  once  again 
to  maintain  the  time  separation  between  P„  and  P,,  the  crustal  thickness  of  the  model  would  have 
to  be  increased.  In  our  model,  crustal  thickness  is  48  km,  which  is  consistent  with  estimates  based 
on  gravity  data  (Hu  et  al.,  1989;  Shi  et  al.,  1989).  For  this  study,  we  also  assume  a  plane-layered 
velocity  structure,  which  would  certainly  be  invalid  if  we  were  to  consider  propagation  from 
earthquakes  farther  south  on  the  Tibetan  Plateau.  Surface-wave  tomography  (Wu,  1993)  suggests 
that  crustal  thickness  varies  only  slightly  in  the  region  from  the  Qaidam  Basin  to  WMQ.  However, 
the  wave  propagation  along  our  profile  may  be  slightly  up-dip,  so  crustal  thickness  should  be 
considered  a  lower  bound,  while  P„  velocity  (8.0  km/s)  is  an  upper  estimate. 

Synthetic  frequency-wavenumber  integration  seismograms  are  plotted  on  the  right  sides  of 
Figures  18  and  19.  The  source  parameters  assumed  are  listed  in  Table  3.  These  include  Harvard 
CMT  mechanisms  (PDE)  when  available;  otherwise  an  average  thrust  mechanism  is  assumed.  Once 
again,  the  source  comer  frequency  has  little  effect,  so  the  synthetics  are  simply  bandpassed  from 
the  ground  displacement  impulse  response  computed  to  4  Hz.  For  these  profiles,  source  depth  is 
assumed  to  be  20  km.  In  Figure  18,  only  P„,  Pg  and  S,  (or  Lg)  are  identifiable,  and  the  synthetics 
provide  a  reasonable  agreement  in  character  to  the  observed  waveforms.  In  Figure  1 9,  the  travel-time 
curves  P,  P„,  PmP,  and  the  surface  reflections  of  P„  and  PmP  are  plotted.  This  is  not  to  suggest  a 
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one-to-one  correspondence  of  arrivals,  but  simply  to  suggest  which  phases  are  predicted  to  arrive 
between  P„  and  P^.  At  the  closest  sources,  P  and  S  are  the  largest  amplitude  arrivals.  The  Moho 
critical  distance  is  about  1 10  km,  so  Pm?  becomes  a  significant  arrival  shortly  beyond  this  range. 
At  the  largest  ranges  (700  km  and  beyond),  the  P„,  wavetrain  consists  of  a  series  of  low-amplitude 
arrivals.  In  the  synthetics,  Pg  begins  with  the  direct  P  wave  travelling  nearly  horizontally  through 
the  crust.  At  the  largest  ranges,  the  observed  Pg  is  more  emergent,  undoubtedly  due  to  scattering 
and  lateral  heterogeneity  in  the  crust. 

Comparisons  of  observed  seismograms  for  each  event  with  synthetics  computed  at  10,  20 
and  30  km  are  shown  in  Figures  21-24.  The  observed  displacements  forevents  88273  and  88320 
(both  239  km.  Figures  21  and  22)  are  very  simple,  with  a  large  arrival  preceded  by  a  small  phase 
that  forms  a  shoulder  on  the  first  peak.  In  our  model,  these  events  are  beyond  critical  range,  so  the 
small  first  arrival  is  P„.  The  large  arrival  is  likely  a  combination  of  the  direct  P  wave,  PqP  and  PmP, 
all  of  which  arrive  close  together  for  a  shallow  source.  In  the  synthetics,  however,  P„  is  distinct 
from  this  larger  arrival,  suggesting  that  the  critical  range  in  our  model  should  be  at  slightly  larger 
range.  The  synthetics  also  predict  a  large  secondary  arrival  corresponding  to  sP^P,  which  moves 
out  with  source  depth.  For  event  88320  (Figure  22),  if  the  arrival  observed  at  about  12-13  sec 
corresponds  to  this  phase,  the  source  depth  must  be  less  than  10  km.  Similarly,  the  arrival  observed 
between  this  phase  and  the  first  arrival  may  be  pP„  or  pPmP.  so  we  interpret  the  source  depth  to  be 
about  7-8  km  (denoted  by  the  arrow;  the  PDE  depth  is  33  km).  Event  88273  (Figure  21)  has  no 
distinct  observed  secondary  arrivals,  but  the  coda  following  the  first  arrival  has  a  duration  of  5-6 
sec.  If  this  coda  consists  of  the  surface-reflected  phases,  the  source  of  this  event  must  be  shallower 
than  10  km  as  well. 

P„  for  event  87356  (317  km.  Figure  23)  is  emergent  and  very  difficult  to  identify.  Perhaps 
appropriately,  the  synthetics  for  this  event  contain  some  wrap-around,  so  P„  is  difficult  to  identify 
in  the  synthetics  as  well.  The  first  substantial  arrivals  predicted  consist  of  a  combination  of  PmP 
and  the  direct  P  wave.  If  we  align  the  synthetics  on  the  predicted  P„  arrival  and  align  these  larger 
arrivals  with  the  first  series  of  large  arrivals  in  the  observed  waveform,  we  may  identify  some  of 
the  later  arrivals.  The  PDE  depth  for  this  event  is  21  km  (denoted  by  the  arrow).  At  this  depth,  the 
second  sequence  of  arrivals  (beginning  at  about  1 8  sec)  is  the  second  crustal  multiple,  (PmP)2-  pPmP 
may  be  interpreted  as  the  relatively  low-frequency  arrival  between  these  two  sequences  of  arrivals. 
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The  first  arrival  from  event  87056  (700  km.  Figure  24)  is  easily  identified  as  P„,  but  the 
observed  waveform  is  much  more  complicated  than  the  synthetics.  This  is  particularly  the  case  for 
the  direct  P  wave,  which  is  large  and  distinct  in  the  synthetics,  but  difficult  to  identify  in  the  observed 
record.  This  is  not  surprizing,  since  the  horizontally  traveling  P  wave  to  a  range  of 700  km  encounters 
substantial  scattering  and  lateral  variations  that  are  not  included  in  the  plane-layered  structure  model 
of  the  synthetics.  Nevertheless,  if  we  identify  the  first  arrival  as  Po,  the  secondary  arrival  at  about 
16  sec  as  P^P  and  the  envelope  of  arrivals  beginning  at  about  22-23  sec  as  the  combination  of  P 
with  the  crustal  multiples,  we  find  that  these  arrival  times  are  consistent  with  the  PDF  depth  of  26 
km  (denoted  by  the  arrow). 

CONCLUSIONS  AND  RECOMMENDATIONS: 

Broad  band  P„-Pg  waveforms  from  eanhquakes  at  regional  distances  recorded  at  station  WMQ 
in  western  China  can  be  well  modeled  using  frequency-wavenumber  integration  synthetics  and 
plane-layered  velocity  structure  models.  Although  not  a  primary  purpose  of  this  study,  we  have 
developed  velocity  structure  models  appropriate  for  regional  wave  propagation  from  three  azimuths 
toward  WMQ.  Common  features  include  a  near-surface  velocity  gradient  (modeled  as  a  lower 
velocity  layer),  a  fairly  constant  velocity  (6.0  -  6.25  km/s)  through  most  of  the  crust,  a  lower-crustal 
layer  or  gradient,  and  a  relatively  thick  crust  (48  -  59  km).  The  model  for  the  profile  of  earthquakes 
SSE  of  WMQ  has  a  thinner  crust  (48  km)  than  that  for  the  NNE  profile  (59  km),  which  is  consistent 
with  the  observation  of  Mangino  and  Ebel  (1992)  that  the  Moho  dips  to  the  NNW  at  WMQ.  On 
the  other  hand,  our  models  (and  others  from  nearby  regions,  such  as  Quin  and  Thurber,  1992)  based 
on  regional  P-wave  modeling  have  a  much  larger  velocity  contrast  at  the  Moho  than  the  teleseismic 
receiver  function  models  of  Mangino  and  Ebel  (1992). 

Broad-band  P„-Pg  waveforms  at  ranges  up  to  1000  km  are  result  from  the  interference  of  a 
variety  of  phases.  Since  some  of  these  depart  upward  from  the  source  and  others  depart  downward, 
this  interference  pattern  is  quite  sensitive  to  source  depth.  This  study  has  shown  that,  in  addition, 
the  interference  pattern  is  also  dependent  on  radiation  pattern,  since  the  important  phases  in  the 
P„-Pg  waveform  are  different  for  different  azimuths  and  different  mechanisms.  For  the  SW  profile, 
waves  that  begin  as  upward-departing  S  waves  (e.g.,  sP„,  sPmP)  dominate  the  waveform.  For  the 
NNE  and  SSE  profiles,  direct  P  and  crustal  multiples  (PmP.  (Pm^Iz)  seem  to  be  most  important. 
What  this  means  is  that  if  accurate  Green’s  functions  can  be  computed  (and  Green’s  functions  are 
only  as  accurate  as  the  velocity  structure  model  assumed),  regional  P„-Pg  waveforms  can  provide 
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excellent  discrimination  for  source  depth  and  mechanism.  In  order  to  be  applied  to  nuclear  event 
discrimination,  it  will  be  necessary  to  "calibrate"  a  region,  by  determining  a  velocity  structure  model 
and  comparing  synthetic  and  observed  waveforms  for  a  source  with  known  depth  and  mechanism. 
Thereafter,  the  Green’s  functions  for  that  path  can  serve  to  discriminate  earthquakes  from  explo¬ 
sions. 
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Table  1  -  Earthquakes  Recogled  at  WMQ 


Date 

Time 

(GMT) 

Lat. 

rN) 

Lon 

CE) 

R 

(km) 

Az 

o 

Depth 

(km) 

mb 

SW  Profile 

10/6/87  (87279) 

1306:20.3 

43.44 

88.55 

82.0 

302 

32 

4.8 

12/17/87  (87351) 

1217:25.0 

41.94 

83.20 

421.9 

59 

33* 

5.1 

8/5/87  (87217) 

1024:21.0 

41.36 

82.11 

534.1 

57 

33* 

4.8 

1/5/87  (87005) 

2252:46.5 

41.96 

81.32 

559.6 

66 

17 

5.9 

1/24/87  (87024a) 

0809:21.0 

41.53 

79.32 

731.2 

67 

29 

5.9 

1/24/87  (87024b) 

1340:40.0 

41.44 

79.25 

740.5 

66 

33* 

5.2 

6/8/87  (87159) 

1330:36.0 

39.79 

14.69 

1175.0 

63 

10 

5.1 

4/30/87  (87120) 

0517:37.0 

39.76 

74.57 

1178.3 

63 

8 

5.7 

NNE  Profile 

10/6/87  (87279) 

1306:20.3 

43.44 

88.55 

82.0 

302 

32 

4.8 

9/18/87  (87261) 

2159:15.0 

47.02 

89.66 

387.6 

204 

33* 

5.3 

4/1/88  (88092) 

0127:16.0 

47.53 

89.64 

439.4 

201 

10 

4.6 

7/23/88  (88205) 

0738:9.7 

48.71 

90.56 

586.4 

203 

19 

5.5 

6/30/88  (88182) 

1525:15.5 

50.23 

91.14 

759.4 

202 

33* 

5.0 

9/16/87  (87259) 

1759:30.6 

52.09 

95.70 

1095.5 

216 

33* 

4.8 

SSE  Profile 

10/6/87  (87279) 

1306:20.3 

43.44 

88.55 

82.0 

302 

32 

4.8 

9/29/88  (88273) 

0700:3.1 

41.75 

88.47 

238.9 

345 

33* 

4.7 

1 1/15/88  (88320) 

1656:46.2 

42.02 

89.30 

239.4 

327 

33* 

5.0 

12/22/87  (87356) 

0016:39.04 

41.36 

89.64 

316.9 

330 

21 

5.9 

2/25/87  (87056) 

1957:52.0 

38.10 

91.18 

699.9 

336 

26 

5.7 

12/6/87  (87340) 

1620:44.9 

37.39 

94.52 

917.8 

323 

33* 

4.7 

Depth  uncertain.  33  km  is  the  PDE  default. 

Compiled  from  PDE,  Wu  (1990),  and  Bennett  et  al.  (1990). 
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lablg  2  z  Yeloci0  Strucmrg  Modgls 


Vp 

(km/s) 

Vs 

(km/s) 

Density 

(g/cm^) 

Thickness 

(km) 

Qp 

Qs 

SW  Profile  (Tienshan  and  Tarim  Basin) 

4.80  2.77  2.58 

9.0 

300 

150 

6.25 

3.61 

2.79 

32.0 

800 

400 

7.25 

4.18 

3.00 

15.0 

1000 

500 

8.00 

4.62 

3.33 

20.0 

1200 

600 

8.10 

4.68 

3.36 

20.0 

1200 

600 

8.20 

4.73 

3.40 

40.0 

1200 

600 

8.30 

4.79 

3.45 

h.s. 

1200 

600 

NNE  Profile  (Altai  Mountains) 

4.80  2.77 

2.58 

5.0 

300 

150 

6.00 

3.46 

2.79 

42.0’ 

800 

400 

6.90 

4.00 

2.85 

45.0^ 

52.0^ 

4.0’ 

1000 

500 

7.40 

4.27 

3.00 

3.0^ 

4.0' 

1.0’ 

1000 

500 

8.15 

4.70 

3.20 

6.0^ 

1.0' 

10.0 

1200 

600 

8.26 

4.77 

3.30 

h.s. 

1200 

600 

SSE  Profile  (near  Lop  Nor) 

4.80  2.77 

2.58 

12.0 

300 

150 

6.25 

3.78 

2.79 

26.0 

1000 

500 

6.80 

3.93 

2.85 

10.0 

1000 

500 

8.00 

4.62 

3.34 

10.0 

1200 

600 

8.10 

4.68 

3.36 

h.s. 

1200 

600 

'  Assuming  a  source  depth  of  10  km. 
^  Assuming  a  source  depth  of  20  km. 
’  Assuming  a  source  depth  of  30  km. 
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Table  3  -  Source 


Date 

R 

(km) 

Az 

o 

Strike 

O 

Dip 

0 

Rake 

D 

Mo 

(xlO^  dyne -cm) 

C 

(Hz) 

SW  Profile 

87279 

82 

302 

220*’ 

40" 

65" 

30." 

Syn2 

200 

60 

220'’ 

40" 

65" 

30." 

87351 

422 

59 

220" 

40" 

65" 

30." 

0.5 

87217,  87005 

560 

66 

lie 

21* 

47* 

41.* 

0.8 

87024a,  87024b 

740 

67 

268* 

45* 

107* 

200.* 

0.8 

Syn6 

950 

60 

220" 

40" 

65" 

30." 

87159 

1175 

63 

298‘ 

27* 

91* 

6.5* 

0.3 

NNE  Profile 

87279 

82 

302 

120" 

90" 

180" 

1.8" 

Syn2 

200 

200 

120" 

90" 

180" 

1.8" 

87261,88091 

400 

200 

154* 

90* 

180* 

1.3* 

88205 

586 

203 

331* 

66* 

170* 

9.0* 

88182 

759 

201 

152* 

50* 

-50* 

3.5* 

Syn6 

950 

200 

70" 

70" 

-40" 

1.8" 

87259 

1095 

216 

70" 

70" 

-40" 

1.8" 

SSE  Profile 

87279 

82 

302 

270" 

60" 

70" 

1.8" 

Syn2 

160 

330 

270" 

60" 

70" 

1.8" 

88273,  88320 

240 

330 

270" 

60" 

70" 

1.8" 

87356 

317 

330 

316* 

53* 

54* 

2.1* 

Syn5 

500 

330 

270" 

60" 

70" 

1.8" 

87056 

700 

336 

16T 

60* 

68* 

5.8* 

87340 

918 

323 

270" 

60" 

70" 

1.8" 

*  Mechanisms  are  Harvard  CMT  solutions  published  in  the  PDE. 

No  mechanism  published.  These  values  are  assumed. 

'  If  not  specified,  is  assumed  to  be  >  4  Hz. 
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Earthquake  Profiles  Recorded  at  WMQ 


1  -  Shaded  topographic  map  of  the  northwestern  border  region  of  China  showing  the  locations 
of  earthquakes  (stars)  recorded  at  CDSN  station  WMQ  (triangle).  Earthquakes  SW.  NNE 
and  SSE  of  WMQ  were  modeled  in  this  study  (see  Table  1).  Also  shown  are  the  locations 
of  the  Kazakh  test  site  and  the  Lop  Nor  test  site  (circles).  Topography  (courtesy  of  Eric 
Fielding,  Cornell)  is  plotted  with  lighter  shades  indicating  higher  elevations.  Gray  squares 
indicate  missing  elevation  data.  The  small  gray  oval  SE  of  WMQ  is  the  Turfan  Basin,  which 
is  below  sea  level.  The  light  region  to  the  south  is  the  edge  of  the  Tibetan  Plateau. 
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Fig.  2  -  Profiles  of  observed  (left)  and  frequency-wavenumber  integration  synthetic  (right)  vertical-component  displacement  waveforms 
for  earthquakes  from  the  SW  recorded  at  CDSN  station  WMQ.  Also  shown  are  travel-time  curves  for  important  P-  and  S-wave 
phases  computed  for  a  source  at  20  km  depth  in  the  velocity  structure  model  listed  in  Table  2. 
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Structure  Model  for  SW  Profile 


Density  (g^cm^)  Velocity  (km/s) 

1  2  3  1  2  3  4  5  6 


80  ^ 


Fig.  4  -  Velocity  structure  model  for  the  profile  SW  of  WMQ. 


87279  (82  km  Range) 


Fig.  5  -  Comparison  of  the  observed  Pa-Pg  waveform  for  event  87279  (top  trace)  with  F-K  synthetics  computed  for  source  depths  of  10, 
20  and  30  km.  All  traces  are  aligned  on  the  P  arrival.  The  moveout  of  significant  phases  (indicated  by  lines  superimposed  on 
the  synthetics)  changes  the  interference  of  arrivals.  This  enables  a  determination  of  .source  depth  (denoted  by  the  arrow  labeled 


ICC)  wilh  F-K  synthetics  computed  for  source  depths  of  10, 
!  traces  arc  aligned  on  the  P,  arrival. 


87005  (559  km  Range) 


Comparison  of  the  observed  P„-P,  waveform  for  event  87005  (top  trace)  with  F-K  synthetics  computed  for  source  depths  of  10. 
20  and  30  km.  The  format  is  the  same  as  in  Figure  6. 


67024a  (731  km  Range) 


Comparison  of  the  observed  P„-P,  waveform  forevent  87024a  (top  trace)  with  F-K  synthetics  computed  for  source  depths  of  10, 
20  and  30  km.  The  format  is  the  same  as  in  Figure  6. 


87159  (1175  km  Range) 
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Fig.  1 1  -  The  portions  of  the  observed  (left)  and  synthetic  (right)  seismogram  profiles  NNE  of  WMQ. 


structure  Models  for  NNE  Profile 


Density  (g/cm^)  Velocity  (km/s) 


Fig.  12  -  Velocity  structure  model  for  the  profile  NNE  of  WMQ. 


87279  (82  km  Range) 


Comparison  of  the  observed  P„-Pg  waveform  for  event  87279  ^top  trace)  with  F-K  synthetics  computed  for  source  depths  of  10, 
20  and  30  km.  The  format  is  the  same  as  in  Figure  5. 


87261  (400  km  Range) 


-  Comparison  of  the  observed  P„*Pg  waveform  for  event  87261  (top  trace)  with  F-K  synthetics  computed  at  400  km  range  for 
source  depths  of  10, 20  and  30  km.  The  format  is  the  same  as  in  Figure  6. 


88092  (400  km  Range) 


-  Comparison  of  the  observed  P„-Pg  viravcform  for  event  88092  (top  trace)  with  F-K  synthetics  computed  at  400  km  range  for 
source  depths  of  10, 20  and  30  km.  The  format  is  the  same  as  in  Figure  6. 
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Fig.  20  -  Velocity  structure  model  for  the  profile  SSE  of  WMQ. 


88273  (239  km  Range) 


Time  (sec) 

Comparison  of  the  observed  P.-P,  waveform  for  event  88273  (top  trace)  with  F-K  synthetics  computed  for 
20  and  30  km.  The  format  is  the  same  as  in  Fivinr 


88320  (239  km  Range) 


Fig.  23  -  Comparison  of  the  observed  Pn'Pg  waveform  for  event  87356  (top  trace)  with  F-K 
20  and  30  km.  The  format  is  the  same  as  in  Figure  6. 
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